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AN OPTIMAL INITIAL GUESS GENERATOR FOR ENTRY 
INTERFACE TARGETERS 

Juan S. Senent* 


If a pure numerical iterative approach is used, targeting entry interface (El) con- 
ditions for nominal and abort return trajectories or for correction maneuvers can 
be computationally expensive. This paper describes an algorithm to obtain an op- 
timal impulsive maneuver that generates a trajectory satisfying a set of El targets: 
inequality constraints on longitude, latitude and azimuth and a fixed flight-path an- 
gle. Most of the calculations require no iterations, making it suitable for real-time 
applications or large trade studies. This algorithm has been used to generate initial 
guesses for abort trajectories during Earth-Moon transfers. 

INTRODUCTION 

Calculating a return trajectory that targets a very specific set of (El) conditions can be computa- 
tionally expensive if a pure numerical iterative approach is considered. The number of constraints 
and the nature of them, as we will see later, can make the problem very difficult for a numerical 
optimizer since multiple local optima and feasibility problems will appear. As an intermediate so- 
lution we can calculate an initial guess that can be used later by the numerical optimizer. This paper 
describes an algorithm to generate such initial guess with the following features: 

• It is analytical, in order to compute the optimal maneuver no numerical iteration is needed. 
The only exception is the calculation of the velocity magnitude from the time of flight. Since 
we are only generating an initial guess, it is beneficial to do it in a very fast way. The nu- 
merical optimizer will obtain the final trajectory taking into account more realistic gravity 
models. 

• It finds the global optimal solution. Depending on the inequality constraints, this optimiza- 
tion problem might present several local optima. This algorithm searches for solutions in all 
the feasible regions guaranteeing a global optimal solution. Algorithms using a numerical 
optimizer might run into a local optimum. 

• It can handle inequality constraints in some of the El targets. Sometimes the El targets are 
defined as a large database of points that create an area. We can check individually each 
point, which can be computationally very expensive, or create one or more rectangular areas 
that include them. In this way we will reduce the computational time dramatically. 

• It can provide information about the feasibility of the problem. One the disadvantages of using 
a pure numerical approach (e.g. a numerical optimizer) is that we do not have an insight of 
the problem. If the optimizer fails to converge, we will not know the reason. It might be due 
to our problem formulation, or an error in the program that runs the simulation or it might be 
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that the problem is not feasible. This algorithm can generate information about the feasibility 
of the problem since it can provide which inequality constraints could not be satisfied. We 
can use the information to redesign the return trajectory or create a new El area. 


We can specify the El targets in different ways: longitude, latitude, azimuth or cross-range and 
down-range, etc. The type of targets and constraints that the algorithm can handle are described in 
Table 1. In the first part of the paper, we will describe the solution to the problem when the inequality 
constraints are not active. In this way, given the state of the spacecraft at a given time, a time of 
flight between the burn and El, and a given velocity magnitude, we will obtain an optimal impulsive 
maneuver that targets: altitude and flight-path angle. We can formulate an equivalent problem by 
using the velocity magnitude instead of the time of flight. Once the basic optimal problem is solved 
then we will add one by one the latitude, longitude and azimuth inequality constraints in the next 
sections of the paper. Adding these inequality constraints will generate disjoint El regions that can 
be reached and therefore each region will have its own local optimum that has to be checked. 


Table 1. El variables and constraints 


El variable 

Description 

Altitude ( h ) 

Fixed 

Longitude (A) 

Inequality constraint 

Latitude ( L ) 

Inequality constraint 

Velocity magnitude (u) 

Fixed (or calculated from time of flight) 

Azimuth ( Az ) 

Inequality constraint 

Flight-path angle (7) 

Fixed 


The applications of this algorithm as an initial guess generator include: a nominal or an abort 
return trajectory from the Moon that targets a specific geographical area at El and that takes into 
account a safe module disposal (avoiding some regions, e.g. some islands) or computing correction 
maneuvers before El to target a very specify El point. This algorithm has been extensively used 
and validated in the generation of abort return trajectories generated after the translunar injection 
maneuver. 


BASIC ALGORITHM 

In this section we will examine the basic El targeting problem with no inequality constraints. The 
statement of the problem is the folowing: 

Given the state of the spacecraft: position and velocity r a , v a , calculate: a minimum impulsive 
maneuver that targets the entry interface position magnitude r e , the entry interface velocity magni- 
tude v e and the entry interface flight-path angle 7 e . Once we solve this problem we can formulate a 
similar one where the target is the time of flight instead of v e . 

Parti 

In order to solve the problem, the algorithm will take several steps, the first part will consist of 
calculating the basic parameters of the return trajectory: 


• Calculate the semi-major axis and eccentricity of the return trajectory (also the parameter of 
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the orbit) 


o _ 


/ i {r e V e COS7 e ) 2 , (r e U e COS7 e ) 

e = \/l p = a(l-e) = 

V ar/ [i 


Calculate the velocity mangnitude after the return manevuer 

w ° = = ^ (h - h) 

where v e is the minimum u e that is necessary to obtain a return trajectory 


Calculate the true anomaly at El 


6 P = arctan 


tan 7 e 


Calculate flight-path angle after the return maneuver 7 a . There are two solutions: 


/ v e v e cos 7e \ 
7 ai ,2 = ± arccos — J 


Additionally, the following condition should be satisfied: 

|r e u e cos7 e | 


Calculate the true anomaly after the return maneuver 0 a for the previous two solutions 


0 ai 2 = arctan 


tail 7ai,2 
1 - ^ 


Calculate the TOF from the abort maneuver to El. First calculate the eccentric anomaly at El: 


E e = ± arccos 


Although there are two solutions to the above equation, the El is before periapse so only the 
negative one will be used. The same equation should be applied to the eccentric anomaly 
after the abort maneuver: 

_ ( a — r a \ 

E ai 2 = ± arccos I ] 


Again, there are solutions to this equation only if 


— <!<=> a(l — e) < r a < a( 1 + e) 


Finally, the TOF can be calculated: 


TOF\2 = U ~ tai'2 = \- Ee ~ Eai ’ 2 ~ e ( sin Ee ~ sin (3) 

Note that the TOF calculation is only valid if e < 1. Therefore this version of the algorithm 
is only applicable for elliptical return trajectories. 
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Two families of solutions have been found. Although these families will satisfy the constraints they 
will have two different TOFs. 


Part II 

In the second part of the algorithm we will calculate the relationship between the position a 
velocity vectors at El: r e , v e . If we use the / and g functions in terms of the transfer angle (A0) 
(see for example ****BOND****) 


r a = r e f + \ e g (4) 

v a = r J + v e g (5) 


where 


/ 

9 

f 

9 


1 (1 — cos AO) 

P 

r a r e . 

— — sin A 6 
s/PP 

— (1 — cos AO) — — 
P r e 

1 — — ( 1 — cos AO) 


HI ( 1 — cos A 0\ 
V p V sin A 6 ) 


1 

r a 


where A 6 = 0 a — 0 e . We can now compute the angle a between the r a and v e : 


r I v e = r J v e/ + v e9 = r e v e sin7 ef + v\g = r a v e cos a 


or alternatively, 


a = arccos 


r e siii 7 e / + v e g 


r a 


r a x V e I = l( r e/ + v e</) x v e | = |/r e x v e | = v e v e cos 7e |/| = r a v e sma 


tana = 


r e v e CQS7e|/| 
r e v e sin 7 e/ + v^g 


r e cos 7e | / 1 
r e sin 7e/ + V e g 


We can see that the angle a is constant and therefore all the possible v e are contained in a cone 
around r a (see Figure 1) whose equation can be defined as: 


v e = Cr a + R cos Au + R sin Av 


( 6 ) 
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Figure 1. Relationship between the r a and v e vectors. 

wherer a is a unit vector in the direction of r a , and u. v will complete the right-handed coordinate 
system. Their definition will change depending on the application. Also the angle A will determine 
the location of v e in the cone {A 6 [0, 2k]). From Figure 1 we obtain (note: a 6 [0, k] ): 


C = v e cos a 
R = v e sin a > 0 


We can define a unit vector in the direction of v e 


v e = cos or a + sin a cos Au + sin a sin A\ 


From eq. 4 we can obtain an expression for r e : 
7 


= J {r a ~ V e £/) 


using eq. 
r, 


*>*> 


= -y — y ( Cr a + R cos Au + R sin Av) 
r ^ C )f 0 -^( C o»ifl + smAv) 


(7) 

( 8 ) 
(9) 
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Again, all the possible values of r e are contained in a cone around r a . 


?e = ( \ r a - ^ (cos Au + sin Av) 

\ J r e / J r e 

Using eq. 9 we can obtain: 



I f r e 

= sin a 


we can express r e in spherical coordinates as (1 don't know if I need this, the second part is much 
better): 


r e = cos pr a + sin p cos Au + sin p sin Av 


where: 


cos p 
sin/? 


r a - gC 
f r e 

_gR 

f r e 


Relationship between p and A^: 


r a — gC \ „ gR , . ^ 

r e = | —7 r a - — (cos Au + sin Av) 

fr e J f r e 


r J r a = cos A 6 — > cos A 6 = cos p 


gR_ 

fr e 


^ sin A9v e sin a 

f r e 

r e cos 7e I/I 


— sin A 6 


xfrpf 


— sin A 9 


m 

/ 
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Therefore: 


sin p = sin A 6 if / < 0 p = A 9 if / < 0 

sin p = — sin A 0 if / > 0 p = — A 9 if / > 0 

if / < 0 — > r e = cos A#r a + sin A 9 cos Au + sin A 9 sin Av 


if / > 0 — > r e = cos A#r a — sin A 9 cos Au — sin A9 sin Av 
If we define f s = — -y-, (f s = ±1) then 


= fs sill A9 

f r e 


and r e can be defined by: 


r e = cos A6r a + f s (sin A 9 cos Au + sin A 9 sin Av) 


and therefore, 


r J v e = cos A 9 cos a + f s sin A 9 sin a cos 2 A + f s sin A 9 sin a sin 2 A 
= cos A# cos a + f s sin A 9 sin a = cos (A 9 — f s a) = sin 7 e 


cos (A# — f s a ) = siii 7 e — > sin ^ — A9 + / s o;^ = siii 7 e 


7e = ~ - + f 3 a 


A 9 


7T /7T \ 

- / s a = 2 - le — ^ ■ sin (A 9 - f s a) = sin - 7 e J = cos7 e > 0, since 7 e G - 


therefore 


sin (A 9 — f s a ) > 0 

Also, if / < 0, we can calculate the angular momentum of the abort transfer orbit as: 


L X V. = 


r a 

cos A 9 
cos a 


u 


f s sin A 9 cos A 
sin a cos A 


v 

f s sin A 9 sin A 
sin a sin A 


( 10 ) 

(ID 

( 12 ) 


7T 7T " 

2 ’ 2 . 

( 13 ) 
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r P x v, 


e = r a ( fs sin A# cos A sin a sin A — sin a cos Af s sin A 0 sin A) 

— u (cos A 0 sin a sin A — cos af s sin A 0 sin A) 

+v (cos A 0 sin a cos A — cos af s sin A 0 cos A) 

= sin A (sin A 0 cos af s — cos A 0 sin a) u — cos A (sin A 0 cos af s — cos A 0 sin a) v 
= f s sin (AO — f s a) (sin ylu — cos Av) 


r e x v e | = sin (A 0 — f s a ) = cos7 e 


h = 


r„ x v. 


sin (A 9 — f s a) 


=f s (sin Au — cos Av) 


(14) 


The angular momentum vector remains in the u — v plane, that is, r Jh = 0 (esto era obvio pues 
r a is on the transfer plane). 

The trajectories can be classified into two families: 


1. The ones that cross the antipode (-r a ) and then the entry interface point (r e ) 

2. The ones that cross the entry interface point (r e ) and then the antipode (-r a ) 


We can find a criteria to distinguish between these two families based on the angular momentum of 
the transfer orbit. We first get — f a x r e 


-fa X r e = 


fa 

-1 

cos A 9 


u 

0 

f s sin A 6 cos A 


v 

0 

f s sin A# sin A 


— r a x r e = — u(— f s sin AO sin A) + v(— f s sin AO cos A) 

= f s sin AO (sin Au — cos Av) 

= sin A0h 

In this way, we can classify the orbits according to this criteria as follows: 

f Family 1 if sin AO > 0 
\ Family 2 if sinA0 < 0 

Also, since 

f a x r e sin AO - - 

77 VV — ~T~- TV 77 h = ±h 

|r a x r e | |sinA0| 

The angular momentum vector is parallel (or antiparallel) to the intersection of the planes defined 
by r a and r e (cosa que ya sabiamos pues los dos vectores pertenecen al piano de transferencia) 


8 



We can now calculate the maximum and minumum inclination of the abort trajectory. If the 
coordinate system is defined: 

11 = |proj- B z| = co^DECr Projection of the z-axis on the plane defined by r a . Where projf a z = z — (: 
v = r a x u = cos *pffe r To complete the right-handed coordinate system 
r a In the direction of r a 

Note: \projt a z\ = cos DEC Ta 

Esto esta MAAAAAALLLLLL. Rehacer con lo que tengo en mis notas: 19 / 01/09 

Since the inclination is the angle between z and h, the minimum inclination will occur when h is 
parallel to u and the maximum inclination of the abort trajectory will occur when h is antiparallel 
to u. We can calculate the angles A associated with these two situations: 

Fromeq. 14 we can get: 


INC min = DEC ra «- A = ~ 

INCmai = 7T — INCmin = 7T — DEC Ta «- A = 

where DEC Va is the declination of r a in the x - y - z coordinate system (intertial system) 
See fig. dibujar figura!!!! 

Some properties that can be obtained: 

(fa X ?e) T V e = fj (f e X V e ) = COS7 e f Jh = 0 


(f a X f e) T projf e v e = (r a X r e ) T v e - (rjv c ) r e = 0 


(r a x r e ) J pro j ia v e = (fa x r e ) T v e - (rjv e ) r a j = 0 

Therefore v e , projt e v e and projt a v e \ie on the same plane (the one defined by r„ x r e ) 

Using eq. 14 we can prove that the amount of plane change only depends on the increment on A. 

If we obtain the abort trajectories defined by A\ and A2, the angle between the two angular 
momentum vectors is: 


hi = sin A 1 11 — cos Aiv 


I12 = sin A2U — cos A2V 


h} I12 = sin A\ sin A2 + cos A\ cos A2 = cos A A = cos 7/ — * A A = 7/ 
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therefore the amount of plane change is defined by the increment in A 
Using eq. 14 we can prove that r 0 and — r G are contained in the plane of the orbit: 


±rjh = ±rj (sin^4u — cosAv) = 0 

We can now calculate v a velocity after the abort maneuver. Using eqs. 15, 7 and 6: 

Va = r e f + v e g = | (r a — v e g) + v e g = ^r a + (g - v e = j (/r a + v e ) 


If \ R 

= 7 r a + C r a + — (cos + sin Av ) 

\ J / / 


(15) 

(16) 


Once again, v a lies on a cone defined by r a 

If v~, velocity vector before the abort maneuver, is known, we can compute the optimal velocity 
after the abort manuverv*. The optimal velocity is just the projection of v^on the cone defined by 
eq. 15. We need to calculate the optimal angle A * that will produce v*. 


V* - L +c\r I R P r °jr« v « 

“ - \f +C ) + f \projr^a 


V a“ V a = T 


proj ra Va _ ( cosylfi + sinylv) 


\P ro Jr a Va 


~ T / * \ 1 1 

11 (v a -v a ) = - 


R 


aT 


\proj ra Va 


R 


v T K - v a ) = - 


aT 


V proj ra v~ _ sinA 
Proj ra V a 


aT, 


cos A* = 


* u proj r a y a 


| proj ra Va 


T, 


sin A* = 

| proj^Va 


tan A* = 


V ProjraV a 
u J proj Ta Wa 


(17) 


Note: if / < 0 then from eq. 15 the radius of the cone j < 0 and therefore the equivalent cone 
will be such as the original but where the uand v are rotated by 7 r. In that case: 


A' = arctan | ^/ ro ^ v ° | + T 
u'proj ra v a 


Note: If the minimum distance between v a and v a (minimum Av ) is obtained when the angle 
A* is such that the unit vector cos.A*u + sin A*v points in the direction of proj Ta y ~ , deviations 
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from A* will increase the distance between v~ and v a (increase in Av). The maximum distance 
(maximum Au) will be found for a v a generated by A* + 7 r, that corresponds to a vector that points 
to —proj ra v~ . Therefore the Av will monotonically increase when A goes from A* to A* T 7 r and 
will monotonically decrease when A goes from .4* + tt to A*. Therefore, if constraints are impose 
on .4, [-Amin 1 Amax] • 


A* 


A* if A* G [Amin, Amax] 

* Amin if A ^ min < &VA max 
k Amax A <Mmax < &V Arnin 


Note: From eq. 17 we can see that if the coordinate system is defined as in page 9 , A* depends 
only on and r a . That is, v e does not affect A*, therefore from eq. 8 we can say that all the 
optimal abort trajectories will be on the same plane, since angular momentum vector h is the same. 
The transfer plane of the optimal abort trajectories will not depend on the v e . 

We can now consider constraints in the latitude at entry interface. 

From fig. la figura de la minima inclinacion, we can express z in the coordinate frame defined in 
9 as: 


z = sin - DEC ra ) u + cos - DEC ra ^j r a = cos DEC Ta t i + sin DEC Ta r a ( 1 8) 

z T r e = sin (DECei) 

z T r e = cos DEC Ta sin AO cos A + cos A# sin DEC ra 


cos A = 


sin DECei — cos AO sin DEC Va 
cos DEC ra sin AO 


(19) 


So if we consider that at entry interface (El) DECei ~ LATei we can calculate the angle A 
that will produce an abort trajectory with a determined LATei- 

Note: for each LATei there are two solutions to eq. 11: ±Alat E i • The positive one will 
generate positive El azimuths (direct return, INC G [0, |]) and the negative one will produce 
negative El azimuths (retrograde return. INC G [f , 7r] )• We can see this if we calculate the INC 
of the abort trajectory: 


cos INC = z T h = sin A cos DEC Ta (20) 

If A G [0. 7r] then we obtain a direct return: 


sin A > 0, since cos DEC Ta > 0 — > cos INC > 0 — » INC G 
If A G [— 7T. 0] then we obtain a retrograde return: 



Azei € [0, 7r] 
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sin A < 0, since cos DEC Ta > 0 — > cos INC < 0 — ► INC € 


TC 

12 '* 


Azei € [tt, 2tt] 


Note: the relationship between the angle A and Azei is given by: 


cos INC = cos LATei sin AZei = sin A cos DEC r 


sin AZei = 


sin A cos DEC ra 
cos LATei 


We can compute AZei as a function of A. 


cos AZei 


(r e x z) 
|r e x z| 


r c x z defines a plane with constant longitude that contains r e 
From their definition we can compute r e x z: 


( 21 ) 


f e xz = — cos DEC Va sin A0cos/ir a +sin DEC Va sin A9 sin >lu+ (cos DEC Ta cos A 0 — sin DEC Ta sin Ad cos ^4) v 
also 


|f e x z| = sin Q - LATei ) 


= cos LATei 


so now we can get 

fr-r (fe X z) 

|f e X Z| 

sin 2 A sin DEC Va sin AO + cos 2 A sin DEC Va sin AO - cos A cos DEC ta cos AO 

cos LATei 

sin DEC Va sin AO — cos A cos DEC Va cos AO 
cos LATei 

Using the above eq. and eq. 21 we can compute tan AZei 


cos AZei = 


tan AZei = 


sin A cos DEC Ta 

sin DEC Ta sin AO — cos A cos DEC Va cos AO 


sin A 

tan DEC Ta sin AO — cos A cos AO 


If a specific AZei is required, we can now solve for A 
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sin A + cos A cos AO tan AZei 


tan AZei tan DEC Ta sin A 6 
tan AZei tan DEC Va sin A 6 


sin (A + e) \/l + cos 2 AO tan 2 AZei — 


where 


, ’ cos AO tan AZei \ 
e = arctan I I 


Therefore we can have two solutions for A 


. . ( tan AZei tan DEC Tn sin AO 

A i = arcsin 

\ V 1 cos 2 AO tan 2 AZei 


— e 


A 2 = 7r — arcsin 


tan AZei tan DEC Ta sin AO 
+ cos 2 AO tan 2 AZei 


— e 


As an alternative: 


sin ( A\ + e) = sin (n — Ai — e) = sin (A 2 + e) 


7r — A\ — e = A 2 + e —> A -2 = 7r — A\ — 2e 
AZei can be also directly computed as follows: 

If we define an E-N-R coordinate system associated with r e : 


R 

N 

E 


fe 

projve Z 
\projr e z 


N x R = 


Z X r e 

\P r °jr e 2 


AZei = arctan 


= arctan 


(jP ro jf e V e ) T ( z X ^e) 
[proji e \ e ) J pro ji e z 
vj (z x r e ) 


= arctan 


(v e -fJv e f e ) T (zxf e ) 


L(v e -fjv e f e ) (z — rjzr e )j 


[z T (v e -fjv e f e )j 



'vj (z x r e )' 


= arctan 

c V / 

“I" 

= arctan 


v^projt e z 



(r e x v e ) 


z projf v e 


Since once the LATei is specified, A is determined using eq. 11, the above eq. is only useful 
to obtain the value of the abort trajectory AZei but not to specify it. That is, we cannot specify 
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LATei and AZei of the abort trajectory at the same time. From the two solutions of the equation 
only one is valid?? averiguar. 

From eq. 12 we can see that the sign of A and AZei&xq the same. 

Calculating the optimal abort trajectory when latitude constraints [LATmin , LAT max ] at El are 
imposed. In the rest of the document we will assume that DECei ~ LATei 

We can rewrite eq. 11: 


cos A = 


1 

cos DEC ra sin A 9 


sin (L AT ei) 


cos A 6 sin DEC Ta 
cos DEC ra sin AO 


= Ki sin (LATei) + K 2 


( 22 ) 


Given the latitude constraints above, a range for the angle A, [ A m i n , A max ], can be calculated: 

First we will calculate the range of feasible latitudes for this problem. The range of feasible 
latitudes [LATfmin, LATf max ]w\\\ be such that that eq. 14 has a solution and also it has to be 
within the range [LAT m in, LATmax ] 


Lmin I\\ sin LAT m in T E 2 


Lmax K 1 sin LATmin + #2 


[j LATfmin j LATf m ax] can be calculated using the following table: 



Lmin > 1 

— 1 ^ Lmin A 1 

T-H 

1 

V 

c 

1 

Lmax ^ 1 

No solution 

[LATmin 1 LATp] 

[LAT m , LAT p ] 

1 A Lmax A 1 

[LAT p , LAI max] 

[LATmin-, LATmax] 

[LATm, LATmax] 

Lmax < — 1 

[LAT P , LAT m ] 

[LATmin-, LAT m ] 

No solution 


where LAT p and LAT m are defined as: 


K\ sin LATei + K 2 = 1 — > LAT P = axcsin ( — — 

V Ki 


I\ \ sin LATei + K 2 = 1 —* LAT m — arcsin ( — — 

V Ki 

Once we have the range of feasible latitudes that the abort trajectories can achieve, we have to 
calculate the range [A m i n , A ma x}' 

Since LATei € [~ f , f ] , hacer grafica de mis notas 7/ 1 4/08 



K\ > 0 — » K\ sin ( LATei)+K 2 increases monotonically 


{ 


A min = arccos (Ki sin (LATfmin) + K 2 ) 
Amax = arccos (Ki sin (LATf max ) + k 2 ) 
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s _ . „ f Amin = arccos (A 'i sin (LATf mc 

K, < 0 - K, sin (LATei)+K 2 decreases monotomcally ^ = azaxa \ Kl sin \ LATfrln 

Note: if the abort trajectory is retrograde then the range would be [— A max , — A m i n ]. 

Algortihm to compute the optimum abort trajectory when latitude constraints are imposed 

1. Specify type of solution 1 or 2 (from eq.3) 

2. Specify direct or posigrade 

3. Compute the optimum A* using eq. 10 

4. If A* € [A m i ni A max \ then END 

5. Compute and AvA max and return the minimum 

We can also find the relationship between A a required longitude at entry interface: 

We can calculate first the difference between the longitude of the antipode of r a and the longitude 
at any entry interface point r e 

(pr°j z r e ) T (proj z - r a ) = \proj z r e \ \proj z - r a | cos AA 

\proj z r e \ = yj 1 - (z T r e ) 2 = 1 - sin 2 ( DECei ) = cos DEC ei 

\proj z - r a | = cos DEC Va 
(proj z r c ) T {proj z - f a ) = -vJproj z t a 


(. proj z r e ) J ( proj z - ? a ) = 


-rJproj z r a 


- f I [r„ - (z T ?a) z] = -fjf a + (z T f e ) (z T f a ) 
— cos A 0 + sili DECei sin DEC Ta 


cos AA = 


— cos A# 4- sin DECei sin DEC r 


cos DECei cos DEC Ta 
So if a specific A A is required, we can compute DECei associated with it: 

sin DECei sin DEC Ta — cos DECei cos DEC Ta cos AA = cos A 6 
sin {DECei + e)^sin 2 DEC ra + cos 2 DEC Va cos 2 AA = cosAfl 


) + K ‘i) 
) + k 2 ) 
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7r — DEC Ei x — t = DECei 2 + e — ► DECei 2 = tt — DEC Eh — 26 

Alternative formulation !!!!!!!!!!!! 

We can also find the relationship between A a required longitude at entry interface: 

We can calculate first the difference between the longitude of the center of the entry interface 
points c s r G and the longitude at any entry interface point r e 

{proj z r e ) T (pro j z c s r a ) = \proj z v e \ \proj z c s t a \ cos A A 
\proj z r e \ = \Jl - (z T r e f = yj 1 - sin 2 ( DEC E i ) = cos DEC E i 
\proj z c s Y a \ = cos DECr a 
(proj z r e ) J ( proj z c s r a ) = c s rJproj z r a 


(proj z r e ) T (pro j z c s r a ) = c s rJproj z r a 

= c s rj r a - (z T r«) zj = c 8 r Jr a - c s (z T r e ) (z T r a ) 

= c s cos A0 — c s sin DECei sin DEC Va = —c s ( — cos A 0 + sin DECei sin DEC Ta ) 


cos AA = 


— c s (— cos A 6 + sin DECei sin DEC Ta ) 
cos DECei cos DEC Ta 


So if a specific AA is required, we can compute DECei associated with it: 
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cos DECei cos DEC ra cos AA + c s sin DECei sin DEC , 


c s cos A 0 


c s cos DECei cos DEC Ta cos A A + sin DECei sin DEC Ta = cos A 6 


sin DECei sin DEC ra + c s cos DECei cos DEC ra cos AA = cos A 6 


sii \{DECei + e 


sin 2 DEC Ta + cos 2 DEC ra cos 2 A A 


cos A 6 


e = arctan 



cos DEC Ta cos A A 
sin DEC Ta 


= arctan 


c s cosAA \ 
tan DEC Ta ) 


Now we can obtain the two solutions for DECei- 


DEC ee = arcsin ^ 

DECei 2 = 7r — arcsin 
sin (DEC Eh +e) = 


cos A^ 


y/ sin’ 2 DEC Ta + cos 2 DEC Ta cos 2 AA 


— e 


cos A 6 


>/sin 2 DEC Va + cos 2 DEC Va cos 2 AA 


— € 


sin (7r - DEC Eh ~ e ) = sin ( DECei 2 + e) 


7r — DEC Eh ~ e = DECei 2 + e —* DECei 2 = 7r — DEC Eh ~ 2 e 

Once we have the DECei associated with the required AA, we can compute the angle A using 
eq. 19. Therefore in order to compute A such that the abort trajectory has a specific longitude at 
entry interface the following procedure should carried out: 

1 . Compute the longitude associate with the antipode of r a 

2. Compute AA = A e — A_ ra 

3. Compute the latitude(s) DECei associated with this value of AA 

4. Compute the angle(s) A assciated with DECei 
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We can also calculate the angle A that can generates the maximum AA. The maximum AA occurs 
when a meridian is tangent to the projection of the circle generated by r Q (see Fig. HACER LA 
FIGURA!!!!!!!): 


(proj^e z) T (proj fa r e ) = 0 


Woj Cs ? a f e 


r e - (c s rjr e ) c s f a ] = [r e - cos A0r a ] 
f s sin A 6 (cos Au+ sin Av) 


z T r e = sin DEC Va cos A 0 + f s sin A 0 cos DEC Ta cos A 


z — ^z T f e ^ r e j [f s sin A 6 (cos Au+ sin Av)] 
|z — ^z T f* e ^ r e (cos Au+ sin Av) 
cos DEC Va cos A — (z T f e ^ (fs sin A 0 cos 2 A + f s sin A 6 sin 2 A) 


cos DEC ra cos A — (sin DEC Ta cos A 0 + f s sin A 0 cos DEC Ta cos A) f 3 sin A 6 
cos DEC ra cos A — sin DEC Ta cos A 0 f s sin A^ — sin 2 A( 9 cos DEC Va cos A 
cos DEC Ta cos A (1 — sin 2 A0) — sin DEC ra cos A 0f s sin A 0 


0 

0 since f s sin A 0 ^ 0 (A^ ^ 

0 

0 

0 

0 


f s sin DEC rn cos A 9 sin A0 

cos J\ — — 

cos DEC Ta (1 — sin 2 A#) 


f s tan DEC ra tan A 6 


If there are no solutions to the eq. it will mean that the maximum AAis reached when A = A m * n 
oi A = A max . 


If a range of longitudes is specified, we can calculate the range(s) of the angle A associated with 
it. 

NOTE ON AZIMUTHS 

In general, we can classify all the cases according to three criteria: family, El plane position and 
orientation of the r e in the (u — v) plane. 

• As it is mentioned above, the family (1 or 2) criteria can obtained by using sin AO 

• El plane position. For the same velocity vector, the azimuth value will change its sign if the 
El plane is on the side of the antipode or not. The criteria we will use is the sign of cos AO. 
If this sign is negative the El plane is on the side of the antipode (see eq. 12) and viceversa. 
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Table 2. Azimuth case classification 


Family 

El plane position 

r e in the (u — v) plane 


Family 

El plane position 



sin A 9 

cos A9 

f s sin A 9 

Case 

sin A 9 

cos A 6 

fs 

Case 

<0 

< 0 

<0 

G 

< 0 

< 0 

>0 

G 

<0 

< 0 

>0 

C 

< 0 

< 0 

<0 

C 

<0 

> 0 

<0 

F 

< 0 

>0 

>0 

F 

<0 

> 0 

>0 

B 

< 0 

>0 

<0 

B 

>0 

< 0 

<0 

E 

> 0 

<0 

<0 

E 

>0 

< 0 

>0 

A 

> 0 

<0 

>0 

A 

>0 

> 0 

<0 

H 

>0 

>0 

<0 

H 

>0 

> 0 

>0 

D 

>0 

>0 

> 0 

D 
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Figure 5. Azimuth type 4 
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The orientation of r e in the (u— v) plane can be easily obtained using f s sin A 6. If f s sin A 6 > 
0 ....EXPLAIN!!!!! 


Some cases are equivalent, for example: cases A and or B and H (see above figure). In some cases 
Az changes continuously when the angle A is varied, in some other cases there is a discontinuity 
when A = 0. Also, the signs of A and Az can be the same or opposite. In this way, 


if fs > 0 — > sign(Az) = sign(A) 

if f s < 0 — * sign(Az) ^ sign(A) (23) 


if f s cos 6 > 0 — ■> discontinuity at A = 0 
if f s cos 6 < 0 — > continuity at A = 0 


Finally, we can classify all of the azimuth types according to the following table: 


fs 

cos A# 

Azimuth type 

>0 

< 0 

1 

Az 

/ 

/ 

A 




<0 

>0 




2 


\^^A 

>0 

> 0 


Az 


3 

A 

A 

<0 

< 0 

4 

Az 



/ 
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Note on the extremes of the arccos function 


y = arccos (/(a)) 


-f'(a 

v/wR 


f"(a) (1 - /(a) 2 ) 1/2 - Hi ~ /H 2 ) V2 (~ 2/(a)/'H) 
1 - f(af 


if y' = 0 <=> f'(a) = 0 , therefore 


/|„'=o = -/"(a)(l-/(a) 2 ) 1/2 


if f"(a) >0 — > y 1 < 0 — > maximum (24) 

if f"(a) <0 — > y" > 0 — > minimum 

Calculating the maximum or minimum azimuth 

From the law of cosines we can get: 

cos l = cos a cos r + sin a sin r cos x 

cos l — cos a cos r 

cos x = : : 

sin a sin r 


where: 

r € [0, tc]: angle between the entry interface position and the antipode 
l G [0, 7r]: angle between the north ( 2 -axis) and the antipode 
a € [0, 7r]: angle between the north ( 2 -axis) and the entry interface position 

From the results of the previous section we can now calculate the maximum or minimum x. We 
only need to calculate the conditions for cos.t. Since r and l are constant, cos a; will depend only 
on a, 

d cos x sin 2 a cos r sin r — cos a sin r (cos l — cos a cos r) cos r — cos a cos l 
da (sin a sin r) 2 sin 2 a sin r 

therefore 


d cos x 
da 


= 0 — » cos r — cos a cos l = 0 — > cos a = 


cos r 
cos l 
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We have a necessary condition for finding the extreme: 


if we define c s 


cos A 9 
| cos A0| 


cos / = z T (c s r a ) = c s sin DEC Ta 
sin l = | z x (c s f a )| = cos DEC ra 
cos r = rj (c s r a ) = |cosA#| 
sinr = |r e x (c s r a )| = |sin A6\ 


tan / = 


cos DEC ra 
c s sin DEC Ta 


{ 


cot DECr a — > l = 5 — DEC ra 
- cot DEC Ta —>/ = | + DEC Va 


if c s > 0 
if c s < 0 


tanr = 


| sin A0| 
cos A0| 


cosr 
cos / 


< 1 


cos A# 
sinD£'C rQ 


< 1 


d 2 cos x 
da 2 


sin a cos / sin 2 a sin r — 2 (cos r — cos a cos l) sin a cos a sin r 
(sin 2 a sin r) 2 


d 2 cos x 
da 2 


d cos x n 

da ’ 


sin 3 a cos / sin r _ sin a cos l 
(sin 2 a sinr) 2 s ^ nr 


If we take into account: 


sinr = |sin A9\ > 0 


sin a > 0 


cos l = c s sin DEC Ta 

From eq. 24 we can obtain a criteria for maximum or minimum azimuth 


d 2 


cosx 


da 2 


^ 0 it c s sin DEC Va ^ 0 * x max 
< 0 it c s sin DE E ra <c 0 > vc-min 


We can find now the maximum or minimum x. First, 
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sin a 


>0^ \/T 


— cos 2 a = 



V cos 2 l — cos 2 r 
|cos/| 


cos X = 


cos l sill r 


cos l — gjy cos r |cos l | (cos 2 l — cos 2 r) |cos l | \/cos 2 / — cos 2 r 

sin a sin r \/cos 2 / — cos 2 r sin r cos / 

| sin DEC ra | \/sin 2 DEC ra — cos 2 A# 
c s sinD£'C ro |sinA#| 

>/ sin 2 DEC Ta — cos 2 A 6 
| sin A# | 


= sign 


, c s sin D.EC ra 


x = arccos 


. ( \/ sin 2 DEC ra - cos 2 A 9 . 

sign ( A/1 , , c s sin DEC Ta 


| sin A9\ 

Finally, Azcan be calculated (see diagrams ....explain) 


f 7 t — x if c s sin A9 < 0 
\ x if c s sin A9 > 0 


if c 3 sin A9 < 0 
if c s sin A 9 > 0 


{ 

{ 


x min 

-> Az max 

X max ~ 

Az m i n 

x min ~ 

Az m i n 

X max ~ 

AZmax 


Also the angle y associated to Az can be calculated as follows: 


cos a = cos l cos r + sin / sin r cos y 


cos y = 


cos a— cos /cos r 
sin l sin r 


— cos / cos r cos r — cos 2 1 cos r 
sin l sin r sin l sin r cos / 

cos r sin/ _ |cos A9\ cos DEC Va _ c s |cosA#| 
sin r cos / |sin A9\ c s sin DEC Va |sin A0| tan DEC Va 


cos A 9 

sin AO | tan DEC ra 


From eq. 23, therefore we can calculate A as: 


y = arccos 


cos A# 


| sin A0| tan DEC Va 


The magnitude and sign of Acan be calculated (see diagrams ....explain ...the sign Acan be deter- 
mined by the signs of Az and f s ) 
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4 _ f sign(7r - y, f s Az) if f s sin AO < 0 
\ sign(y , f s Az ) if / s sin AO > 0 

The above derivations calculate only Az £ [0. tt] if we are interested in retrograde solutions 
(Az ' , A'), all we have to do is: 


Az' = -Az 
A' = -A 


Az m in * Az max 
Az max —> Az min 


CALCULATING V E GIVEN THE TOF 


e cos E e 
e sin E e 



a 


r e v e _ r e i> e sin7 e 

y/JId y/gd 


Note: 


= ^y m 

y/gd y/gd 


ar e v e 


- yari = = y / Te (2a - r«) 

y/gd 


= _ i = ?p(2 ^ ' P ~ _ x = ( 2 « ~ r p) _ 1 = Q~r p 


r e u e sin 7 e _ r e v e sin 7 e _ siii7 e \/r e (2a - r e ) 
e> /ua a-r p 


COsE e = 


a — r e 
ea 


a - r e 


a-r p 

a 


a 


a — r e 
a — r p 


tan E e = 


r e Ve sin 7 e 
y/JIa 


1 - 


Te_ 

a 


areVe siii7e . / — r 

yfiia _ sin 7e y r e (2 a - r e ) 
a — r e a — r e 


cos E a 


a — r 


a 


ea 


a ~ r a 
a — r p 


L = r e — a 
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M = — sin 7 e y/ r e (2 a - r e ) 


I\ = r a — a 


sin E e = 


-M 
a — r p 


cos E e 




L 2 

1 — sin 2 E e 


L 2 

M 2 

(a— r p ) 2 




(a - r p ) 2 L 2 
(a - r p ) 2 - i \/ 2 


1 = 


L 2 


(a — r p ) 2 — M 2 


(a - r p y 


2 rn _ ( a ~ r a\ 2 _ K 2 

ms (^ a _ rp J - /,2 + M 2 

I ^ 2 “ 

““ = T-ZHF 

The sign of sin E a will depend on the type of solution: 1 or 2 

if E a G (-7T, 0) Ea < 0 — ► sin E a < OTypel 


if E a G (0. 7r) E’a > 0 — ► sin E a > 0Type2 


cos E a = 


a — r a 


-K 


a -r p sjL 2 + M 2 


Note: 


a — r p > 0 


that's why we get the positive solution of \/L 2 + M 2 


tan E a 


± yj l ~ u+m 

-K 

y/L 2 +M 2 


± VL 2 +M 2 -K 2 
-K 


-r p ?L 2 
r p f - M 2 

= L 2 + M 2 
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